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ABSTRACT 


An algorithm used to composite SSM/I 85.5 GHz imagery and derive sea ice 
motion was adapted for operational testing at Fleet Numerical Meteorology and 
Oceanography Command (FNMOC). A feature tracking technique was applied to a 6- 
month period, with data provided by FNMOC and the Naval Research Laboratory (NRL). 
Ice motions are detectable using the SSM/I motion algorithm, and fields of SSM/I motion 
vectors are qualitatively consistent with coincident fields of in situ buoy motion vectors. 
Accuracy of the SSM/I motion vectors relative to buoy motion vectors increase 
significantly with buoy speed. No correlation between SSM/I and buoy motion vectors 1s 
observed for speeds below 3 cm/s and correlation increases significantly above 5 cm/s. 
The results are very sensitive to compositing techniques used to combine SSM/I passes 
into a single sea ice representation. FNMOC data was composited using a “drop-in-the- 
bucket" technique while NRL data was composited by a bi-linear interpolation technique. 


Significantly poorer results were found with FNMOC composited data. 
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I. INTRODUCTION 


Passive microwave observations of polar oceans have become essential to the 
tracking of ice motions, estimating sea ice concentrations, and for classifying sea ice 
types. Ice formation and movement has always been a concern to Naval, DOD and 
commercial shipping. Global ice information, immediately practical for use in shipping 
and petroleum development activities, has broader implications for the meteorological 


foundations of modeling climate change. 


Remotely sensed microwave imagery provides an excellent data source for 
computation of sea ice motions. Resulting ice motions can further be assimilated into 
polar ice prediction models used by the National Ice Center (NIC) and Fleet Numerical 


Meteorology and Oceanography Center (FNMOC). 


Buoys set in Arctic ice provide a continuous, although spatially limited source of 
ice motion information. Buoy positions are accurate to within 350 meters (Colony and 
Thorndike 1984). Though sparse and unevenly distributed, they are useful in evaluating 
the quality of SSM/I derived ice motions. The scarcity of in-situ buoy observations leads 


to the importance of using satellite-based SSM/I analysis to derive ice motion. 














A. MOTIVATION 


Previous studies of derived SSM/I ice motion have centered on using archived 
data. Temporal periods have ranged from 24 and 72-hour derived ice motions, to 
monthly and multi-year climatic studies (Emery et al. 1991, Emery et al., 1997, Kwok et 
al., 1998). Traditionally, ice motions are derived for the entire Arctic basin. Computing 
daily SSM/I ice motions will provide FNMOC and NIC real-time ice motions that are 


useful for application in daily ice motion products and various modeling programs. 


Knowledge of the location and extent of sea ice has always been of vital 
importance. Assimilating sea ice motions into polar ice prediction models would 
certainly produce better operational ice motion products. Our only direct method of 
obtaining surface meteorological parameters in the Arctic comes from a limited number 
of buoys spread throughout the Arctic Ocean. They provide highly accurate temporal and 
spatial data, though very sparse in spatial coverage. Ice motions produced on a synoptic 
scale will provide useful input for future ice motion product development at NIC and 
modeling development for the U.S. Navy's Polar Ice Prediction System (PIPS 3.0) model. 
Incorporating daily SSM/I sea ice motion data should greatly enhance accuracies of these 


products. 


Ice motions in this study are derived from passive microwave imagery acquired 
by the Special Sensor Microwave/Imager (SSM/I) aboard the Defense Meteorological 
Satellite Program (DMSP) satellites F-11/13/14/15 (F-11 was turned off Spring 2000). 
The microwave spectrum is used as clouds are mostly transparent and ice has a higher 


emittance than water, making it detectable using the Microwave Imager. Data acquired 











by the Advanced Very High Resolution Radiometer (AVHRR) has a sack higher 
resolution, but visible and infrared radiation is not transparent to clouds. Data acquired by 
the Synthetic Aperture Radar (SAR) is limited by spatial coverage and very large files 
that cause time acquisition constraints that prevent the day-to-day operational use 
required by FNMOC and NIC. The DMSP SSM/I program is designed to provide a 


continuous flow of data, which FNMOC assimilation techniques require for their models. 


B. OBJECTIVES 


The object of this paper is to perform a quantitative and qualitative assessment of 
ice motions derived from an ice motion algorithm developed by Meier (1999). Case 
studies are developed to analyze SSM/I derived ice motions and compare them to in-situ 
Arctic buoy motions. Normalized SSM/I and buoy speed differences, and relative 
differences in SSM/I and buoy direction are computed on a daily timeframe. The 
correlation of buoy and SSM/I motion vectors provide a quantitative evaluation of the ice 
motion algorithm and the value of SSM/I derived vectors. To be useful for this purpose, 
SSM/I derived motions will be analyzed daily. Our goal is to imitate FNMOC’s daily 
routine of receiving the raw SSM/I data, processing it through the ice motion algorithm, 
and making derived ice motions available for assimilation into polar prediction models. 
By imitating FNMOC procedures, we can provide a quality assessment that is useful to 


determine any algorithm characteristics that may require updating for future application. 


Resulting correlations will provide FNMOC and NIC a tool from which to assess 
SSM/I derived ice motion accuracy and aid in assignment of a quality parameter, used for 


assimilating parameters into polar models. 
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Ul. THEORY 


A. MEASUREMENTS 


Satellite multi-channel passive microwave sensors provide global radiance 
measurements from which to map, monitor and study Arctic sea ice motion. The SSM/I 
microwave radiometer measures the emitted energy from the earth and atmosphere in the . 
microwave frequency spectrum. The seven frequencies (in GHz) and polarizations of the 
SSM/I are 19.3 V (vertical polarization), 19.3 H (horizontal polarization), 22.235 V, 37.0 
V, 37.0 H, 85.5 V, and 85.5 H. For sea ice motions, the 85.5 GHz horizontally polarized 
channel is used, since the spatial resolution (12.5 km) of this channel and its sensitivity to 
ice properties provides the best opportunity to compute ice motions. 


B. RADIATIVE TRANSFER 


Radiative transfer theory provides the basis for methods used to determine sea-ice 
motion from remote sensing. In the microwave part of the EM spectrum we typically 
work with frequency (v) rather than wavelength (A) as the spectral variable. The 85.5 
GHz frequency is located in a region of high transmittance relative to lower surrounding 
water vapor and oxygen transmittance bands. Figure 1 shows an annual mean 
transmittance value of 0.8 at 85.5 GHz for 15° north latitude. Transmittance is typically 


greater at polar latitudes due to reduced water vapor within the atmosphere. 
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Figure 1. Frequency/Wavelength Spectrum. 


Frequencies are small enough (or wavelengths long enough) so that the Planck 
function can be simplified into a linear function referred to as Rayleigh-Jeans 
Approximation (Kidder and Vonder Haar, 1995). Temperatures encountered on the earth 
and atmosphere at centimeter and millimeter wavelengths can be estimated. Since the 
Rayleigh-Jeans Approximation holds in the microwave wavelength regime, the emitted 


energy is proportional to the temperature of the radiator. 


The emission characteristics of the passive microwave signature of sea ice are a 
function of electromagnetic properties of the surface and surface temperature. The 
brightness temperature when transmittance is equal to 1.0 is defined as Tg = eT, where € 
is the emissivity and T, is the physical surface temperature. Emissivity is the fraction of 
radiant energy that a perfect black body would emit at the same physical temperature, 
dependent upon the electromagnetic (EM) properties of the emitting surface, frequency 


and polarization. 


The emitted microwave radiance is sensitive to the emission characteristics of 
open water, first-year ice and multi-year ice. Water is highly reflective in the microwave 


range and characterized by strong differences between horizontal and vertical 
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polarization. Emission values of water for the 85.5GHz channel are 0.79 for vertical 
polarization and 0.53 for horizontal polarization (Meier 1998) and therefore appear 
radiatively cold. First year ice has an emissivity value of around 0.9. The higher value 
makes the ice appear radiatively warmer than open water and therefore detectable in 
contrast to water. The emissive value of multi-year ice is variable as properties of old ice 
are complex and its structure varies. The typical value is found to be around 0.7 and 
therefore appears radiatively cooler than first year ice. Complications in interpreting 
microwave signatures occur during the spring and summer months as the 85.5 GHz 
channel is susceptible to atmospheric water vapor and emissive properties of melt ponds 
that collect on the ice surface. This gives the ice a lower emissive value, therefore 


appearing radiatively cold as in open water. 


C, THE 85.5 GHZ SENSOR 


A particular advantage of using microwave frequencies is the insensitivity to solar 
irradiance, clouds and other atmospheric effects that inhibit visible and infrared 
frequencies. The 85.5 GHz microwave frequency sensor is used for ice motion detection 
as it provides better resolution required in determining small ice motions normally found 
in the Arctic. Table 1 provides SSM/I sensor characteristics for all channels. 


Table 1. SSM/I channel characteristics 




























Calculation of ice motions using microwave imagery does not require detection 


and detailed analysis of ice characteristics. The only requirement is that ice features 
dominate the variations of brightness temperature in a given image. Detection of the 
movement of these features in successive images can then be confidently ascribed to ice 


motion. 














Hr. DATA 


To validate the SSM/I ice-motion algorithm for the real time application 
described in this thesis, data sets were selected based on the availability of in-situ buoy 
measurements and receipt of near real time satellite data. Composited satellite images 
with resulting ice motion vectors were matched to buoy observations both in space and 
time. This chapter will describe the data sets and instrumentation used to collect the data 
used in this study. 

A. INSTRUMENTS 

i: Special Sensor Microwave/Imager (SSM/I) 

The SSM/I is a seven-channel, four-frequency, linearly polarized, passive 
microwave radiometric system that measures scene microwave brightness temperatures at 
19.35, 22.235, 37.0, and 85.5 GHz (Table 1.). SSM/I is carried aboard Defense 
Meteorological Satellite Program (DMSP) satellites, F11, F13, Fl4 and F15. They 
maintain an altitude of approximately 830km in a polar orbit. DMSP satellites provide 
excellent polar coverage, as the orbital period is approximately 102 minutes. This is 
especially beneficial for microwave imaging of sea ice, as numerous passes will provide a 


dataset that covers the polar region in a little over 12 hours. 
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Figure 2. Orbital characteristics of the DMSP Satellite and 85.5 GHz channel. 
Figure 2 shows swath data characteristics of the SSM/I imager (PODAAC, 2000). 
A swath consists of A/B scan pairs, which include 256 radiance measurements. Large 
circles signify all SSM/I microwave channels, while the small circles signify the 85.5 
GHz channel. The 85.5 GHz channel provides double the information than other 
channels and at better resolution. A constant 45.4° angle of view measures the upwelling 


scene brightness temperatures allowing wind effects to be detectable on the surface. 


Zi Arctic Buoys 

The International Arctic buoy Programme (IABP) is located in the Polar Science 
Center, Applied Physics Laboratory, at the University of Washington. They maintain a 
network of automated data buoys in the Arctic Basin that monitor synoptic-scale fields of 


pressure, temperature, and ice motion that provide real-time support to Arctic operations 
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and meteorological and oceanographic research. Table 2 provides a listing of buoys 
tracked by the IABP. A total of 28 different buoys were used for this study, most of 
which were ICEX-AIR buoys. Figure 3 provides an example of the ICEX-AIR Arctic 
buoy that provided Meteorological and motion parameters used in this study. Detailed 


buoy diagrams are located on the IABP web page (IABP 1996). 
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Figure 3. An ICEX-AIR buoy capsule. 
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DATE ARGOS WMO 
DEPLOYED ID ID 
Aug 1999 1222 48601 
Aug 1996 1261 48102 
Apr 1998 1301 48581 
Aug 1998 1351 48532 
Aug 1998 1793 48533 
Aug 1998 1905 25525 
Aug 1998 1906 25526 
Apr 2000 

. Sep 1998 2388 25557 
Mar 2000 2415 48579 
Apr 1998 2416 47523 
Sep 1997 2417 48572 
Mar 2000 2418 48577 
Aug 1998 3004 25535 
Aug 1999 3311 25569 
Aug 1999 3691 25012 
Apr 1998 4954 48580 
Apr 2000 5300 
Aug 1999 5312 48523 
Mar 1999 5313 47538 
Mar 2000 5317 48524 
Mar 2000 5318 48525 
Apr 2000 7091 26509 
Apr 2000 7092 26510 
Apr 2000 8059 71541 
Aug 1999 8065 25570 
Apr 2000 8066 71547 
May 1992 10667 48531 
Aug 1996 19578 47602 
Apr 2000 20726 695 
Aug 1999 20858 48574 
Aug 1999 20859 48575 
Apr 2000 22203 
Sep 1997 22204 
Apr 2000 22206 
Apr 2000 22207 
Sep 1997 26696 48576 | 
Jul 1996 26699 48573 





Table 2. Buoy characteristics (May 1", 2000) 


EXPR 
NUMBER 
484 
484 
1053 
484 
29 
557 
557 
557 
1053 
1053 
1053 
1053 
1053 
1053 
919 
314 


1053 
627 
627 
627 
627 
627 
282 
282 
919 
919 
919 
1016 
1053 


1053 
1053 


1053 
1053 


GTS 

HEADER 
SSVX01-LFPW 
SSVX01-LFPW 
SSVX02-CWEG 
SSVX01-LFPW 
SSVX01-LFPW 
SSVX12-KARS 
SSVX12-KARS 
SSVX12-KARS 
SSVX12-KARS 
SSVX02-CWEG 
SSVX02-CWEG 
SSVX02-CWEG 
SSVX02-CWEG 
SSVX12-KARS 
SSVX01-LFPW 
SSVX01-LFPW 


SSVX02-CWEG 
SSVX02-CWEG 
SSVX02-CWEG 
SSVX02-CWEG 
SSVX02-CWEG 
SSVX02-CWEG 


-SSVX16-KARS 


SSVX16-KARS 
SSVX01-LFPW 
SSVX01-LFPW 
SSVX01-LFPW 
SSVX02-CWEG 
SSVX12-KARS 


SSVX12-KARS 
SSVX12-KARS 
695 

695 

695 

695 
SSVX12-KARS 
SSVX02-CWEG 
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Latitude 


78.247 
76.823 
83.496 
79.474 
76.287 
78.106 
80.566 
89.486 
78.754 
79.691 
75.084 
79.202 
77.646 
81.29 

82.444 
77,241 


75.622 
89.58 
75.518 
78.94 
84.803 
81.47 
89.351 
88 
88.804 
77.015 
87.35 
71.672 
4:125 
89.49 
76.062 
77.009 
89.497 
80.146 
89.488 
89.522 
79.606 
84.381 


Longitude 


-178.728 
-177.884 
-79.899 
-144.666 
-150.107 
147.28 
172.59 
107.518 
-157,.273 
-131.255 
-138.695 
-146.843 
-133.176 
154.186 
162.32 
175.92 


-139.509 
-107.301 
-176.566 
-120.64 
-116.723 
~111.34 
-12.784 
161.688 
-91.473 
163.65 
-92.609 
170.519 
-126.329 
-105.661 
-176.6 
-166.904 
-105.525 
~154.58 
-105.332 
-107.314 
-156.326 
-96.966 


Data 
Bytes 
17 

17 

16 

17 

17 
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Buoy 
Description 
ICEX-AIR 
ICEX-AIR 
Met ocean 
ICEX-AIR 
ICEX-AIR 
ICEX-AIR 
ICEX-AIR 
JAMSTEC 
CES 

CES 
ZENO-3200 
CES 

CES 
ICEX-AIR 
ICEX-AIR 
ICEX-AIR 


ZENO-3200 
CALIB 
ICEX-AIR 
ICEX-AIR 
EC 

EC 

NAVO 
NAVO 
AWI 
ICEX-AIR 
AWI 
IOEB 
ICEX-AIR 
PMEL 
ICEX-AIR 
ICEX-AIR 
PMEL 
PMEL 
PMEL 
PMEL 
CES 

CES 











B. DATA SETS 

1. FNMOC SSM/I Data 

Fleet Numerical provided daily SSM/I files from February 15" through May 30", 
2000. A minimum number of case studies were selected for analysis and comparison. 
Exceptionally large SSM/I data files and composited imagery prohibited the storage of 6- 
months worth of data. A subjective process was used to isolate periods of Meteorological 


conditions that the author deemed important for application of Meier’s algonthm. 


7: Naval Research Laboratory (NRL) SSM/I Data 

The premise of this paper was to study a 6-month period of ice-motions. FNMOC 
data only covered February — May, so NRL provided SSM/I data in TeraScan data format 
from December 1*, 1999 through February 29", 2000. This provided an overlap where 
both FNMOC and NRL data were available for 15-29 February. NRL and FNMOC 
apply different geolocation and antenna correction patterns to their SSM/I data sets. 
Therefore, the differences in ice-motion results derived from the two data sets were 
examined during the overlap period. Also, compositing NRL SSM/I data was 
accomplished using a routine in TeraScan. A February case study is examined to 


determine if compositing differences occur and result in different SSM/I motion vectors. 


3. International Arctic Buoy Programme (IABP) 
Arctic buoy data was made available on a monthly basis from IABP. Table 3 


depicts typical buoy data characteristics. This study uses buoy U and V wind 
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components, and latitude and longitude coordinates for purposes of plotting and deriving 


buoy motion vectors. 


Table 3. Example data file from [ABP 


Buo 
ae Year  Julianday GMT Latitude Longitude SLP(mb) Air-Temp© U (cm/s) V (cm/s) 
14955 1999 335 0 85.8052 -169.878 1017 -18.3 0.8287 -5.4236 
14955 1999 §=335 12 85.8028 -169.388 1017 -17.6 1.2372 -11.224 
14955 1999 336 0 85.7984 -168.681 1017 -16.6 1.2094 -9.9548 
14955 1999 9336 12 85.7971 -168.325 1019 -16.7 1.0648 -6.9224 


4. NOGAPS 

The Naval Operational Global Atmospheric Prediction System (NOGAPS) 
Meteorological fields provided a data set of surface temperature, sea-level pressure in 
millibars (mbs) and 10-meter winds at one-degree resolution. NOGAPS fields were used 
as a ground truth to check motions provided by buoy observations and compare 
directional SSM/I vectors against expected motions resulting from a given synoptic wind 


pattern. 
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IV. PROCEDURE 


AG OVERVIEW 


SSM/I ice motion fields we expected to be most accurate during winter, when 
surface and atmospheric conditions are optimum for tracking. The SSM/I derived motion 
fields are most valuable for representing large-scale, persistent patterns of sea ice 
transport. Consecutive 85.5 GHz SSM/I imagery provide comprehensive estimates of sea 
ice motion on a uniform grid for sea ice covered areas of the Arctic. The scale resolution 
is not adequate for comparison of fine-scale motions over shorter time intervals, but is 
suited for investigating coarser-scale patterns and processes affected by synoptic-scale 
weather events. 

B. METHODOLOGY 

1. Synoptic Situation 

Long-term climatic studies of ice motion fields using ice-tracking algorithms 
applied to the 85.5 GHz microwave imagery have produced results representative of the 
general Arctic circulation pattern (Emery et al., 1997). This thesis used daily 14-hour 
composites to derive ice motions on a time-scale reflective of FNMOC's daily operating 
procedures. Meteorological synoptic mean sea level pressure, 10-meter winds, and buoy 
locations were plotted simultaneously with the ice motion analysis. The cases cover 
periods of different meteorological conditions throughout the six-month period. A list of 


Case studies is provided in Table 4. 
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Table 4. Selected case studies 


Case | December 27 - 29" 1999 | Low Pressure transits western Arctic 
with High Pressure and increasing 
pressure gradient. 


2000 





















Low Pressure transits western Arctic 
with distinctive wind shift north of 
Nome AK and increasing pressure 

gradient. 

Light pressure gradient as Arctic High 

Pressure ridge extends from west to 

east. 


Intensifying Arctic High Pressure with 
a strong constant pressure gradient. 














— 


February 17 - 18" 2000 
March 06 - 08" 2000 


Ze Ice Motion Algorithm 
















The ice motion algorithm for 85.5 GHz SSM/I imagery presented here was 
developed by Chuck Fowler and Bill Emery at the university of Colorado (Fowler 1995; 
Emery et al. 1997) and adapted by Meier (1999) for integration at FNMOC. Deriving sea 
ice motions from a 24hr, 72hr or monthly satellite composite may mask the algorithm's 
ability to depict short-term changes. The 14-hour period used here was selected to ensure 
complete geographical satellite coverage of the Arctic Ocean and provide the highest 
possible time resolution. 

‘The microwave imager can detect change in the varying ice characteristics by 
detecting the change in brightness temperature (Tg) features in a target area. The basic 
approach involves comparing two spatially consecutive composite images 12 or 24 hours 


apart and searching for the similar Tg patterns in each image (Meier 1999). The 
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displacement of a Tp pattern between images is determined and a motion vector is 


calculated. 


A contiguous matrix of 10x10 Tg target areas would produce an ice displacement 
calculation every Tg pixel (12.5km) and plots one motion vector every 10 Tg pixels 
(125km). A routine within the algorithm allows the user to overlap a portion of the Tp 
matrix. An overlap of 5 Tg pixels was chosen, resulting in ice motion vectors being 
plotted every 5 Tp pixels (62.5 km). 

Another routine of the ice motion algorithm allows over sampling the array to 
calculate motion. Over sampling moves the target area in 1/4-pixel increments. This 
allows detection of Tg features to be correlated and ice motion to be derived at a 
displacement resolution of 3.1 km or speed resolution of 3.62 cm/s. 

3. Ice Motions 

Four steps are required to produce ice motions from FNMOC SSM/I raw data 
files. 

a. Routines within the algorithm are used to convert raw SSM/I files 
obtained from FNMOC into brightness temperatures. This is a two-step process where 
the data is first converted to antenna temperature (T,) and then converted into brightness 


temperatures (Tp). 
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vectors plotted every 62.5krn. 


Brightness Temperature (Tp) pixels. 
gridded onto the NSDIC polar stereographic grid (NSDIC 1997). 
Language (IDL) procedures are used for compositing and mapping. The 85.5 GHz | 
horizontal polarization channel is extracted from the Tg, files, composited and gridded. 


The IDL sub-procedure mapil.pro was modified as part of this analysis to correct the 


ee 


Figure 4(a-d). 4a is a composited 12.5 km resolution satellite image containing 
608x896 Tz pixels. 4b is a standard ice motion vector plot at every 10 Tz pixels. 
A 5 Tz pixel overlap in 4c computes ice motion vectors every 62.5 km. 4d 
reflects over sampling a Tg pixel at 1/4 increment, computing ice motion every 


3.1 km or 3.62 cm/s. 


b. Figure 4(a) reflects a composite that consists of a 608 x 896 grid of 


gridding routines. 
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Brightness temperatures are composited and 


Interactive Data 





As earlier stated, a 14-hour composite is used for this study. The 
algorithm bins Tg values using the drop-in-the-bucket method (Tp pixel values are binned 
from individual SSM/I passes and averaged). A 6-hour, 12-hour and 24-hour composite 
was investigated. A 6-hour composite left several open gaps in coverage. A 12-hour 
composite often left a small section of the Arctic with incomplete coverage. A 24-hour 
composite increases temporal averaging, which could lead to smoothing and excessive 
binning of brightness temperatures, which skew actual motion results. The 14-hour 
composite provided complete coverage of the Arctic and was chosen to reflect FNMOC's 
12-hour watch routine and reduce motion uncertainties that result from smoothing. 

C Once composite images are created, a motion field is calculated. 
The algorithm performs an automated method of calculating ice motion. A feature match 
between images is determined by the maximum correlation (minimum difference) of Tp 
pixel values in a 10x10 pixel target area. The maximum cross correlation (MCC) of the 
target area provides an objective method for computing ice motion from consecutive 
composited images. Correlations of a target area in the first image are compared to all 
regions in a pre-set search area in the second image. 

Figure 4(b) is a standard ice motion vector plot at every 10 Tg pixels. The 
highest correlated peak (MCC method) within the target search area is located and an ice 
motion derived. One ice motion vector is plotted every 125 km or 14.42 cm/s. Figure 
4(c) shows the 5 Tg pixel overlap used for this study. This allows the computation of ice 
motion vectors every 62.5 km or 7.23 cm/s. Figure 4(d) shows the ‘effect of over 
sampling. The target area is sampled in 1/4 Tg pixel increments, computing ice motion at 


a spatial displacement of 3.1 km or speed resolution of 3.62 cm/s. 
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Filtering of motion vectors is applied in order to delete erroneous motions. 


A pre-set correlation threshold value of 0.7 in the algorithm’s filter routine was initially 
used for this study. Correlation is a measure of the match between consecutive target 
area arrays. Any matching feature that produced an ice motion vector in a target area 
between images that correlated below 0.7 will be filtered and not plotted. 

In theory, a correlation threshold value of 0.9 has a higher risk of filtering 
good motion vectors, just as a 0.7 correlation threshold risks inclusion of a higher number 
of bad vectors. This thesis examines the effect of the correlation filter. Motion vector 
fields were generated with values of 0.7 and 0.9 and the results will be discussed in 
Chapter V. 

d. The filtered motion file is plotted using an IDL routine. Figure 5 
shows an example of plotted motions using this routine. To conduct a comparison 
analysis between FNMOC and NRL SSM/I data with in-situ buoy motions, the ASCII 


file of motion vectors were plotted with TeraScan. 
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Figure 5. SSM/I 85.5GHz motion vector image using IDL. 

Z. SSM/I/Buoy Vector Matching 

Two methods were used in this study to compare SSM/I motion vectors with buoy 
motion vectors. First, the closest SSM/I vector within a specified Tg pixel grid radius 
was selected and matched to a co-located buoy. Figure 6 shows the closest SSM/I vector 
in a target area, matched to a buoy. . In theory, one SSM/I vector will always be located 
within a 2.5 Tp pixel radius (31.25 km) centered on the buoy. The correlation threshold 
of 0.7 or 0.9 will remove many ice motion vectors, resulting in SSM/I vector motion 


gaps. This often leaves a buoy unmatched with an SSM/I motion vector. 
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Selecting the closest SSM/I motion vector does not always provide an 
adequate qualitative assessment of the actual SSM/I motion vector field. Correlated 
motion vectors are often erroneous or sparse, leaving no match for a buoy. The second 
method expands the Tg pixel radius to include more than one SSM/I ice motion vector. 
Figure 7 shows a 5 Tg pixel radius or 62.5 km. Extending our radius enables up to 5 
SSM/I ice motion vectors to be located near, and centered on a buoy. The median value 
of the SSM/I motion vectors is then matched to its corresponding buoy. As with method 


one, the SSM/I ice motion vectors must meet the assigned correlation threshold. 


Radius = 2.5 To pools pr 31.25 km 





Figure 6. Closest SSM/I/Buoy match grid cell dimensions. 
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Figure 7. Median SSM/I/Buoy match grid cell dimensions. 
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V. RESULTS AND DISCUSSION 


A. CASE STUDY DISCUSSION 


The following case studies display full coverage SSM/I and buoy motion vector 
fields. The effects of varying the correlation threshold between 0.7 and 0.9, and 
differences between matching a buoy to its nearest SSM/I motion vector and to the 
median of neighboring SSM/I vectors, are described. 

1. Case Study 1: December 27 - 29", 1999 (SSM/I from NRL) 

Figure 8 is a full coverage plot of SSM/I and buoy motion vectors for the period 
of December 27 - 29", 1999 at 1200Z. A definite synoptic change is reflected as surface 
winds increase in response to a strengthening pressure gradient. SSM/I and buoy motion 
vector magnitudes substantially increase in response to the increased wind speeds. . 

Figure 9 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
correlation threshold for December 27", 1999 at 1200Z. The left-hand plot matches the 
closest SSM/I motion vector, with the right-hand plot matching the median SSM/I motion 
vector to its corresponding buoy. SSM/I and buoy motion vector pairs increase as a 
result of using the median SSM/I motion vector. Overall, SSM/I speed and direction 
values are better reflective of buoy values using the median SSM/I vector. 

Figure 10 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
and 0.9 correlation threshold for December 27", 1999 at 1200Z. Overall, SSM/I speed 


and direction values with 0.9 correlation threshold better reflect buoy values. One good 
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SSM/TI and buoy match using 0.7 correlation threshold was eliminated with the 0.9 
correlation threshold. 

Figure 11 is a SSM/I and buoy motion vector match comparison plot for the 
closest and median matching techniques using a 0.7 correlation threshold for December 
28", 1999 at 1200Z. The closest SSM/I motion vector match does better on this day. The 
median SSM/I vector set creates more pairs and SSM/I speed and direction values are 
distinctly more reflective of buoy values. 

Figure 12 is a plot of SSM/I and buoy motion vector pairs comparing a 0.7 and 
0.9 correlation threshold for December 23h 1999 at 1200Z. Results are similar to 
December 27" in that SSM/I motion vectors with a 0.9 correlation filter agree with buoy 
vector values better than a 0.7 correlation. Two good SSM/I and buoy vector pairs using 
0.7 correlation threshold are eliminated with the 0.9 correlation threshold. 

Figure 13 is a SSM/I and buoy motion vector field image for December 28", 1999 
at 1200Z. The closest neighboring SSM/I motion vector in the left-hand image is 
selected and highlighted by ‘Zoom area 1'. Tg grid cells are visibly detectable on the 
right-hand image and reveal that the only SSM/I vector within a 2.5 Tb pixels radius 
(31.25 km) from its matching buoy has an erroneous vector speed of 29 cm/s and an 
erroneous direction. Lower in the image, the buoy with a motion of 7 cm/s appears to be 
surrounded by accurate SSM/I motion vectors. In this example, the closest SSM/I to 
buoy match approach fails to include information from good vectors surrounding the 
buoy (see Fig. 6). 

Referring back to Figure 11, the median SSM/I motion vector method depicts 
several new SSM/I and buoy motion vector pairs. Expanding the Tg pixel radius to 5 
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(62.5km) and computing an SSM/I median, allows the opportunity to include good 
vectors in our analysis (see Fig. 7). 

Figure 14 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
correlation threshold for December 29", 1999 at 1200Z. The synoptic pattern and high 
surface winds are clearly reflected in the SSM/I and buoy motion vectors. The median 
SSM/I vector match clearly adds value, as more pairs are displayed and earlier SSM/I 
speed and direction values are more reflective of buoy values. Two matched pairs are 
clearly erroneous, but they passed the 0.7 correlation threshold filter and are plotted. 

Figure 15 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
and 0.9 correlation threshold for December 29", 1999 at 1200Z. The 0.9 correlation 
threshold removes one of the erroneous pairs and results in more agreement between 
SSM/I and buoy values for two pairs of vectors. Two good SSM/I and buoy pairs are 
removed. 

Summary for Case study 1 (NRL SSM/I data): The December 27 - 29" 
time period was selected to identify differences in resulting SSM/I motion vectors by 
applying different methods of SSM/I motion vector selections, and by varying the 
correlation threshold. This time period starts with buoy motions of 3-6 cm/s that increase 
to 10-26 cm/s in response to an increasing surface pressure gradient. 

Case 1 provided insightful information concerning the SSM/I motion 
vector selection process. Selection of the closest SSM/I motion vector to the 
corresponding buoy resulted in a small number of buoy pairs, and often paired the buoy 
with an erroneous SSM/I vector. Applying a median to SSM/I motion vectors within a 5 


Tp radius provided better results. Several more SSM/I and buoy motion vector pairs 
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were selected and the qualitative comparisons suggest improvement in their agreement. 


Correlation thresholds of 0.7 and 0.9 mirrored one another as both provided excellent 
SSM/I and buoy motion vector pairs. The 0.9 correlations tended to remove good and 
bad vector pairs, but overall provided slightly better values when compared to the 0.7 


correlation threshold. 
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Figure 8. Plotted SSM/I and buoy motion vectors for December 27, 28 and 29", 1999 at 1200Z. 
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Figure 10. Median SSM/I motion vector match using a 0.7 and 0.9 correlat 
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Figure 11. Closest and Median SSM/I mot 
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Figure 13. Highlighted area of closest SSM/I and buoy motion vector plot 
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Figure 14. Closest and median SSM/I motion vec 
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Figure 15. Median SSM/I motion vector match us 
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2. Case Study 2: January 15 - 16th, 2000 (SSM/I from NRL) 

Figure 16 is a full coverage plot of SSM/I and buoy motion vectors for the period 
of January 15th, 2000 at 1200Z and January 16", 2000 at 0000Z. A low-pressure system 
is clearly detected by SSM/I and buoy motion vectors west of the Queen Elizabeth 
Islands. Very high SSM/I and buoy motion vectors cover an extensive region of the 
Western Arctic as the ice responds to very strong surface winds forced by a strong 
pressure gradient. 

Figure 17 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
correlation threshold for January 15", 2000 at 1200Z. Very few closest SSM/I and buoy 
motion vector pairs result in this case. The median SSM/I motion vector set provides 
several more matched pairs. Cyclonic motion vectors reflect the surface low pressure by 
the Queen Elizabeth Islands. One erroneous SSM/I vector from matching only the 
closest vector is retained using the median matching of SSM/I and buoy motion vectors. 

Figure 18 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
and 0.9 correlation threshold for January 15", 2000 at 1200Z. The 0.9 correlation 
threshold retains most of the good SSM/I and buoy pairs, including the cyclonic vectors 
associated with the low-pressure area. However, several good pairs from the 0.7 
threshold do not pass the 0.9 threshold. 

Figure 19 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
correlation threshold for January 16", 2000 at 0000Z. The closest SSM/I motion vector 
match does very well for this 0000Z run. No significant differences in synoptic 


conditions are evident compared to the 1200Z composite on January 15". More pairs 
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result from the median SSM/I and buoy motion vector matching, including two obvious 
erroneous pairs. 

Figure 20 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
and 0.9 correlation threshold for January 16", 2000 at 0000Z. The 0.9 correlation 
threshold removes all clearly erroneous SSM/I and buoy matched pairs and retains the 
good matches of the remaining pairs. 

Figure 21 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
correlation threshold for January 16", 2000 at 1200Z. The closest SSM/I motion vector 
match does very well. In this case, more pairs result from the median SSM/I vector, 
including some clearly erroneous pairs. 

Figure 22 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
and 0.9 correlation threshold for January 16", 2000 at 1200Z. The 0.9 correlation 
threshold removes all clearly erroneous SSM/I and buoy matched pairs and two 
reasonable matched pairs. Unlike the previous day, where motion values improved with 
the 0.9 correlation threshold, no improvement is apparent in motion values that pass the 


0.9 correlation threshold. 


Summary for Case study 2 (NRL SSM/1 data): The January 15 - 16" time frame 
displays high motions resulting from a constant strong pressure gradient. A closed low 
pressure remains stationary northwest of the Queen Elizabeth Islands. SSM/I motion 
vectors reflect both meteorological conditions very well. SSM/I and Baey motion pairs 
are numerous using all three methods. The composite for January 15" at 1200Z did not 
provide as many SSM/I and buoy motion vector pairs as January 16” at 0000Z and 


1200Z. Further investigation revealed fewer buoys were available for January 15". 
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Figure 17. Closest and median SSM/I motion vector match using a 0.7 correlat 
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3. Case Study 3: February 17 - 18th, 2000 (SSM/I Data from 
NRL and FNMOC) 

Case 3 begins a transition from using solely NRL SSM/I data to FNMOC SSM/T 
data. The most important difference between the data sets is the compositing technique. 
NRL formatted data is composited using TeraScan while FNMOC formatted data is 
composited using the Meier ice-motion algorithm. TeraScan uses a more sophisticated 
bi-linear interpolation technique, while FNMOC uses a “drop-in-the-bucket” averaging 
technique. The same SSM/I orbital passes were selected and processed for both 
composite routines in Case Study 3. 

The synoptic condition for Case Study 3 is vastly different from Cases 1 and 2. 
Buoy motions reflect light winds as the western Arctic is under the influence of a high- 
pressure ridge. Case Study 3 shows algorithm characteristics during a light motion 
regime. A comparison between NRL and FNMOC derived motion vectors are also 
described. 

Figure 23 is a full coverage plot of SSM/I and buoy motion vectors for the period 
of February 18, 2000 at 0000Z. The NRL plot is on the left hand side with the FNMOC 
plot on the right hand side. Buoy motions are very small due to light winds. With buoy 
motions less than 3 cm/s, the accuracy of derived SSM/I motion vectors is not expected 
to be high. 

Differences between the NRL and FNMOC SSM/I motion vector plots are 
apparent in Figure 23. Though buoy motion vectors are small and hard to see due to slow 


motions the NRL vectors show buoy directions well. The FNMOC SSM/I motion 
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vectors are erratic and do not provide a sense of ice movement or agree with actual buoy 
direction. 

Figure 24 is a full coverage plot of SSM/I and buoy motion vectors for the period 
of February 18'", 2000 at 1200Z. The NRL plot is on the left hand side with the FNMOC 
plot on the right hand side. Buoy motions remain low throughout this period. The only 
significant difference noted from 12 hours earlier is fewer SSM/I motion vectors are 
plotted for NRL and FNMOC. NRL SSM/I motion vectors continue to reflect the slow 
buoy motion whereas FNMOC SSM/I motion vectors do not reflect the buoy motion. 

Figure 25 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
correlation threshold for February 18th, 2000 at O000Z. The closest SSM/I motion vector 
pairs are plotted for NRL and FNMOC. No reliable pairs are evident. NRL SSM/ 
motion vectors do reflect buoy direction whereas FNMOC does not. An interesting point 
is that the SSM/I and buoy motion vector pairs for both plots are different. Since the 
same data was used in the compositing routines, the correlated pairs show that the 
different compositing routines produce significantly different motion vectors. 

Figure 26 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
correlation threshold for February 18th, 2000 at 0000Z. The median SSM/I motion 
vector pairs are plotted for NRL and FNMOC. No discernible improvements are 
reflected in the NRL plot. The median matching of FNMOC SSM/I and buoy motion 
vectors results in several new pairs and the values match the buoys better than the closest 
buoy comparison. 

Figure 27 is a SSM/I and buoy motion vector match comparison plot using a 0.9 
correlation threshold for February 18th, 2000 at 0000Z. The SSM/I motion vector pairs 
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are plotted for NRL and FNMOC. No discernible improvements are reflected in either 
plot. 

Summary for Case study 3 (NRL and FNMOC SSM/I data): The period 
of February 17 - 18" was selected to investigate the ice motion algorithm ability to detect 
ice motion during a light wind period. This period af time also provided an opportunity to 
directly compare NRL and FNMOC derived SSM/I motion vector results. As expected, 
SSM/I motion vector speeds are not accurate and directions are erratic due to the light ice 
motion conditions. There are clear SSM/I motion vector differences between NRL and 
FNMOC. Buoy motions are hard to see, but the NRL SSM/I vectors reflect buoy 
directional vectors. FNMOC SSM/I motion vectors are sparse and erratic, with no vector 


component matching with the buoys. 
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4, Case Study 4: March 06 - 08th, 2000 (SSM/I from FNMOC) 

Case Study 4 takes a direct look at composited FNMOC data using the ice motion 
algorithm. FNMOC SSM/I motion vectors in Case Study 3 resulted in no accurate SSM/I 
motion vectors. Case Study 4 is influenced by an intensifying high-pressure system over 
the Western Arctic and increasing surface winds that force high buoy motion speeds. 

Figure 28 is a full coverage plot of SSM/I and buoy motion vectors for the period 
of March 06 - 08th, 2000 at 1200Z. SSM/I motion vectors provide a good indication of 
ice movement in the area of high buoy motions. There are several gaps in coverage as in 
Case Study 3. 

Figure 29 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
correlation threshold for March 07th, 2000 at 1200Z. In contrast to Case Study 3, the 
closest SSM/I to buoy motion vector match does relatively well in areas of faster ice 
motion. SSM/I motion vector speeds and directions show a definite improvement due to 
high buoy motion speeds. The median vector match clearly adds value resulting in more 
pairs and SSM/I speed and direction values are in agreement with buoy values. 

Figure 30 is a SSM/I and buoy motion vector match comparison plot using a 0.7 
and 0.9 correlation threshold for March 07th, 2000 at 1200Z. Most of the SSM/I and 
buoy motion vector pairs are removed using the 0.9 correlation threshold. Only one of 
several apparently accurate matched pairs remain. This is very different from the NRL 
0.9 correlation threshold results. NRL 0.9 correlated SSM/I and buoy motion vector pairs 
showed slightly better values compared to NRL 0.7 correlated plots. The 0.9 correlation 
threshold removes many good pairs from the FNMOC SSM/I and buoy motion vector 
pairs. 
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Summary for Case study 4 (NRL and FNMOC SSM/I data): The period 
of March 06 - 08" had an Arctic high-pressure that intensified to 1060 millibars. A 
section of the western Arctic had small buoy motions, whereas north and west of 
Wrangel Island displayed very high buoy motions. There are large gaps in SSM/I motion 
vector coverage as was also seen in Case Study 3. 

The FNMOC median 0.7 correlation threshold plots provide much better 
coverage in SSM/I and buoy motion vector pairs the median 0.9 correlation threshold. 
The 0.9 correlation threshold removed most SSM/I and buoy motion vector pairs. This is 
not true for the NRL composited data. In Case Studies 1 and 2, the 0.7 and 0.9 
correlations provided good SSM/I and buoy motion vector pairs, with the 0.9 correlation 
threshold providing slightly better values. This suggests that under the FNMOC 


compositing, very few vectors have a correlation greater than about 0.7. 
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B. SUMMARY STATISTICS 


Tables 5 through 8 provide normalized summary statistics derived from the 0.7 
and 0.9 correlation thresholds as applied to the median SSM/I and buoy motion vector 
pairs for December through May. A Normalized Speed (NS) value is defined by dividing 
the SSM/I speed by the buoy speed. A NS value equal to 1.0 indicates perfect agreement 
with the buoy speed. The NS is multiplied by 100% for these tables. The Relative 
Direction (RD) is defined as the relative difference between the SSM/I and buoy derived 


directions. 
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Table 5 contains summary statistics for NRL composited data using the 0.7 
correlation threshold. NS values between 80 - 120% are defined as good. As buoy speed 
increased, so do good NS percentages. For Table 5, there were 447 SSM/I and buoy 
pairs with ice motion speeds between 7.5 - 15.0 cm/s. 34.2% of those pairs fall within 
the 80 - 120 % NS range. Of the total pairs, 40.9 % have normalized speed values greater 
than 200%. Almost all of these bad NS values had buoy speeds between 0 - 7.5 cm/s. 
RD values of +/- 30° are defined as good. As buoy speeds increased, so did the good RD 
percentages. 


Table 5. Normalized Speed and Relative Direction statistics for the 0.7 correlation 
threshold December 01, 1999 to February 29, 2000. (NRL data) 


Total Pairs Buoy Speed (emis) 


Normalized 
Speed 


Relative 
Direction 
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Table 6 contains summary statistics for NRL composited data using the 0.9 
correlation threshold. As buoy speed increased, so do the good NS percentages. For 
Table 6, there were 338 SSM/I and buoy pairs with ice motion speeds between 7.5 - 15.0 
cm/s. 36.4% of those pairs fall within the 80 - 120 % NS range. Of the total pairs, 39.6 
% have normalized speed values greater than 200%. As in table 5, most bad NS values 
had buoy speeds between 0 - 7.5 cm/s. As buoy speeds increased, so did the good RD 
percentage values. 


Table 6. Normalized Speed and Relative Direction statistics for the 0.9 correlation 


threshold December 01, 1999 to February 29, 2000. (NRL data) 


0%) 
0.59 


0 
~ Relative. ~ 
Direction 





As shown in the case study analysis, applying the 0.9 correlation threshold to 
NRL composited data resulted in slightly better SSM/I and buoy motion vector values. 
This is also evident in Tables 5 and 6. Care must be taken when viewing statistics of 


buoy speeds greater than 22.5 cm/s due to the limited number of occurrences. 
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Table 7 contains summary statistics for FNMOC composited data using the 0.7 
correlation threshold. NS values between 80 - 120% are defined as good. NS percentage 
values only increased for the 7.5 - 15.0 cm/s range of buoy speed. There were 221 SSM/I 
and buoy pairs with ice motion speeds between 7.5 - 15.0 cm/s and 26.2% of these pairs 
fall within the 80 - 120 % NS range. Of the total pairs, 45.1 % have normalized speed 
values greater than 200%. RD values of +/- 30° are defined as good. As buoy speeds 
increased, the RD percentages did not always increase for the FNMOC data. 


Table 7. Normalized Speed and Relative Direction statistics for the 0.7 correlation 
threshold December 01, 1999 to February 29, 2000. (FNMOC data). 


> 200% 
Normalized 
Direction 


5(29.4%) 
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Table 8 contains summary statistics for FNMOC composited data using the 0.9 
correlation threshold. NS percentage values for the 7.5 - 15.0 cm/s range of buoy speed 
are worse than those reflected in Table 7. For Table 8, there were 164 SSM/I and buoy 
pairs with ice motion speeds between 7.5 - 15.0 cm/s. Of the total pairs, 39.7 % have 
normalized speed values greater than 200%. Almost all of these values had buoy speeds 
between 0 - 7.5 cm/s. This is 3.6 % worse than using the 0.7 correlation threshold. As in 
Table 7, most of these values had buoy speeds between 0 - 7.5 cm/s. As buoy speeds 
increased, so did good RD percentages. 


Table 8. Normalized Speed and Relative Direction statistics for the 0.9 correlation 


threshold December 01, 1999 to February 29, 2000. (FNMOC data) 


rota Pairs Buoy Speed (emis 































Normalized 0-7.5 7.5 - 15.0 15.0 - 22.5 2225 
547.7%) 
> 200% 4(10.0%) 






Relative. 
Direction 
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There is a large disparity between NRL and FNMOC statistics when applying the 


0.9 correlation threshold. NRL data in Table 5 showed higher percentages than FNMOC. 
A key statistic is the total number of SSM/I and buoy vector pairs that have a NS value 
between 80 - 120 %. NRL had a 19.6% NS percentage value between 80 — 120% and 


FNMOC had a 14.0% NS percentage value between 80 — 120%. The difference in NS 
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percentages result from compositing differences applied to NRL and FNMOC data. As 
shown in case study analysis, applying the 0.9 correlation threshold to FNMOC 
composited data resulted in poor SSM/I and buoy motion vector values. This is also 
evident in Tables 7 and 8. Poor results are likely attributable to the "drop-in-the-bucket" 


compositing technique. 


Figures 31-34 are polar plots used to visualize normalized SSM/I and buoy 
motion vectors. The points on the figure indicate the speed and direction of the SSM/I 
vector relative to the buoy vector. Normalized speed is computed by dividing SSM/I 
derived speed by buoy speed, with 1.00 indicating buoy speed equal to SSM/I speed. 
Relative direction was computed by taking the difference between SSM/I and buoy 
derived directions. A SSM/I vector at 0° and 1.0 normalized speed is equal to the buoy 


vector. 


Scaling is standard logarithmic, centered at .01 and ranging to 100. Gray 
background shading is applied to provide the reader a sense of directional bias and 
highlight regions of good (+/- 30°), questionable (+/- 30° - 90°) and bad (+/- 90° - 180°) 
normalized direction. Buoy speeds are grouped by color every 5 cm/s starting at 0.00 


cm/s, with a +/- 2.5 cm/s range. 
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Figure 31 is a Normalized Polar Plot using a 0.7 correlation threshold applied to 


NRL SSM/I data from December 01‘, 1999 to February 29", 2000. A clear progression 


towards a normalized value of speed = 1 and direction = 0° is reflected as buoy speed 


increases. 
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Figure 31. Polar plot using a 0.7 correlation threshold (NRL data) December 01, 
1999 to February 29, 2000. 
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Figure 32 is a Normalized Polar Plot using a 0.9 correlation threshold applied to 
NRL SSM/I data from December 01%, 1999 to February 29", 2000. Progression to 
normalization as buoy speed increases is clearly depicted with fewer reports plotted in the 


bad direction range, indicating the removal of bad vector pairs. 
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Figure 32. Polar plot using a 0.9 correlation threshold (NRL data) 
December 01, 1999 to February 29, 2000. 
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Figure 33 is a Normalized Polar Plot using a 0.7 correlation threshold applied to 
FNMOC SSM/I data from February 15" to May 30", 2000. A lot of noise is reflected in 
the plot. A clear progression towards normalization as buoy speed increases is not as 
evident. Several high buoy speeds are plotted outside the defined good speed and 


direction range. 
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Figure 33. Polar plot using a 0.7 correlation threshold on (FNMOC data) 
February 15 to May 30, 2000. 
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Figure 34 is a Normalized Polar Plot using a 0.9 correlation threshold applied to 
FNMOC SSM/I data from February 15" to May 30", 2000. Clear removal of several 
SSM/I and buoy pairs occur both outside and inside defined good ranges. The plot is 


very noisy with no clear progression to normalization as buoy speed increases. 
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Figure 34. Polar plot using a 0.9 correlation threshold on (FNMOC data) 
February 15 to May 30, 2000. 
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C. | NRL/FNMOC Statistical summary for May 11-14" 

Case Study 3 for February provided a comparison study between NRL and 
FNMOC data sets. Unfortunately, low ice motions during that time period prevented a 
complete analysis for comparing compositing techniques and the effect on the accuracy 
of ice motions. Ice motions were derived and compared for a three-day period in May 
using NRL and FNMOC SSM/I data. This time frame was chosen because high ice 
motions prevail during a period of strong synoptic-scale surface winds. 

Tables 9 and 10 contain summary statistics applying a 0.7 correlation threshold 
for FNMOC and NRL composited data for May 11-14". The total percentage for good 
pairs (80-120%) of NS is much higher for the NRL data set. NRL ice motions also show 
a strong tendency for the SSM/I vectors to match the buoy vectors as buoy speed 
increases. The fraction of good NS values in the FNMOC data set are significantly lower 
than in the NRL data set. NRL also produces a greater fraction of Relative Direction 


values in the good range (+/- 0-30°) than does the FNMOC data set. 
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Table 9. Normalized Speed and Relative Direction statistics for the 0.7 correlation 
threshold, May 11 - 14‘", 2000 (FNMOC data). 


Total Pairs Buoy Speed (cms) 


Normalized 0-7.5 7.5 - 15.0 15.0 - 22.5 2225 
Speed 114 32 45 20 8 





Table 10. Normalized Speed and Relative Direction statistics for the 0.7 correlation 


threshold, May 11 - 14", 2000 (NRL data). 


7 Total Pars 


Normalized 
Speed 


~ Relative 
Direction 


) 
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Figures 35 and 36 are Normalized Polar Plots using a 0.7 correlation threshold 
applied to FNMOC and NRL SSM/I data from May 11-14". A clear progression towards 


a normalized value of speed equal to 1 and direction equal to 0° is apparent in the NRL 


plot. FNMOC normalized results are poor in comparison to NRL. 
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Figure 35. NRL Polar plot using a 0.7 correlation threshold 
May 11 - 14"", 2000. 
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Figure 36. FNMOC Polar plot using a 0.7 correlation threshold 
May 11 - 14", 2000. 
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D. VECTOR CORRELATION 


A vector correlation technique developed by Crosby et al. (1993) was used in this 
thesis to quantify the quality of SSM/I ice motion vectors. This technique addresses the 
problem of correlating vector quantities, unique to meteorological and oceanographic 
studies. Vectors are correlated directly, without decomposition into scalar components. 
The result is an objective measure of the correlation of two sets of coincident vector 
pairs. 

Code for the Crosby technique was used to compute the square of the vector 
correlation coefficient (R’). An R’ value of 2.0 indicates perfect correlation between the 
vector pairs. Crosby (1993) points out that the sample size must be large enough to 
include sufficient variation among the vector pairs. Therefore, since the value of R? 
depends slightly on sample size, all calculations below are done with a constant number 
of vector pairs. To include all information provided by matched pairs and provide a 
statistically meaningful result, all buoy speeds were collected and sorted starting at 0.00 
cm/s. Groups of 25 SSM/I and buoy vector pairs were sequentially entered into Crosby's 


vector correlation program and plotted in Figures 37 and 38 by their average buoy speed. 
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Figure 37 is a vector correlation plot of NRL and FNMOC SSM/I data, applied to 
SSM/I and buoy motion vector results for the 0.7 correlation threshold. The NRL plot 
reflects a strong increase in R’ beginning around 5 cm/s. The FNMOC vectors have 


consistently low R’ values at all buoy speeds. 
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Figure 37. Vector correlation coefficient (R’) plot of 0.7 correlated threshold. 
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Figure 38 is a vector correlation plot of NRL and FNMOC SSM/I data, applied to 
the 0.9 correlation threshold, SSM/I and buoy motion vector results. The 0.9 correlation 
threshold clearly reflects a transition of the NRL data similar to that seen in Figure 43. 
The 0.9 correlation threshold plot for FNMOC is worse than seen in Figure 43 indicating 


largely uncorrelated vectors. 
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Figure 38. Vector correlation coefficient (R’) plot of 0.9 correlated threshold. 


Two conclusions can be made. First, the 0.9 correlation threshold using NRL data 
provided better results as compared to FNMOC. Most FNMOC data is removed with the 
0.9 correlation threshold, leaving a small number of SSM/I and buoy motion vectors to 


correlate. Second, the oversampling scheme of the ice motion algorithm improves the 
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resolution of the vector speed calculation. Without over sampling, SSM/I speeds aren't 
detectable below 7.23 cm/s and Figure 38 clearly indicates correlated values below 7.23 


cm/s. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


Meier's adapted 85.5 GHz SSM/I sea ice motion algorithm was applied to SSM/I 
data received from NRL and FNMOC. Four case studies were subjectively chosen to 
represent different meteorological conditions and conduct comparisons between NRL and 
FNMOC. Within the algorithms filtering subroutine, 0.7 and 0.9 correlation thresholds 
are applied to the SSM/I motion vector field. To assess the quality of the satellite 
analysis, SSM/I motion vectors were matched with in-situ buoys using two methods: (1) 
the closest SSM/I vector paired with the buoy motion vector and (2) the median of SSM/I 


vectors within 5 pixels paired with the buoy motion vector. 


Initial study of the algorithm revealed a mapping error contained within the IDL 
mapper program code. This program converts the geo-reference and mapping Fortran 
routines into IDL code for adaptation by FNMOC. The code was updated and provided 


to FNMOC. 


As applied to NRL SSM/1 data, the 0.7 correlation threshold and median SSM/I to 
buoy motion vector pairs provided good results. NRL motion vectors correlate well for 
speeds greater than 4-5 cm/s and correlation increases with higher speeds. A vector 
correlation coefficient (R*) was computed using a method developed by Crosby et al., 
(1993). Resulting plots indicate a definitive increase in correlation beginning at about 5 


cm/s. The 0.9 correlation threshold does provide slightly increased vector correlation and 


qd 

















a clearer increase in correlation, though a few good SSM/I and buoy motion vector pairs 


are removed. 


When applying the vector correlation to FNMOC SSM/I data, the 0.7 correlation 
threshold and median SSM/I to buoy motion vector pairs display poor results. The 
Statistical analysis reflects noise and ambiguity, with no distinct increase in correlation as 
in the NRL vector correlation plot. Vector correlation values are much lower than that of 
NRL SSM/I data. The 0.9 correlation threshold provided worse results as most good 
SSM/I and buoy vector pairs were filtered, resulting in yet lower R? values. This 
suggests that the FNMOC composting produced few motion vectors with pattern 


correlation greater than 0.7. 


Compositing methods appear to be the primary cause of the high statistical 
differences between NRL and FNMOC SSM/I data. The compositing scheme applied to 
FNMOC SSM/I data resulted in artificial texture that introduced noise into the motion 
calculation and produced poor vector correlation with in situ buoys. Compositing NRL 
SSM/I data in TeraScan produced better ice motion vectors and improved statistical 
results. 


B. RECOMMENDATIONS 


The results presented here strongly suggest that the FNMOC compositing use a 
bi-linear interpolation technique (as used in TeraScan) in the compositing procedure. 
This will provide the best possible ice motions to the models. A 24-hour or 72-hour 
motion composite will produce smoother results and display more SSM/I and buoy 
motion vectors. But care must be given to masking the algorithm's weakness to handle 


Tg changes noted during the melt-by-day, freeze-by-night springtime conditions. 


78 























The removal of SSM/I vectors greater than 90° from wind direction should be 
considered. Kwok et al., (1998) filters S SM/I derived vectors greater than 90° contrary to 
the gradient wind flow. This can be easily applied, and makes sense when using | to 3 
day composites. Caution must be exercised as ice motion can lag in response to wind 


flow or be influenced by a strong wind gradient 100's of kilometers away. 


Case studies for this paper were selective and subjective. A real-time extended 
comparison study using the current algorithm compositing technique and TeraScan 
composite routine with 0.7 and 0.9-correlation filter for 14/24/72 hour periods should be 
conducted. Application to a full Arctic seasonal cycle (October - May) using NRL and 


FNMOC data may provide better insight to the algorithms strengths and weaknesses. 
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